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ABSTRACT 

In this paper, we study the dynamical stabiHty and time evolution of the central dark matter 
cores of low-mass (about 10** — IO^Mq) galactic haloes found in recent cold dark matter 
simulations at high redshift. From those simulations we extract three haloes, assembled by 
hierarchical merging, that at redshift z ^ 10 display a core and we evolve them without fur- 
ther merging to low redshift using direct A^-body integration. The central core in the dark 
matter profile is found to be dynamically stable: it survives for many crossing times without 
evolution into a cusp. This result supports the claim that the mass dependence of the central 
dark matter profiles of simulated haloes is a direct consequence of the power spectrum of 
primordial density fluctuations. In addition, we show that the simulated dark matter profiles, 
if they evolved in isolation, are consistent with the observed velocity dispersion profile of 
stars in the inner parts of the Draco dwarf spheroidal galaxy. Simple scaling arguments are 
reviewed which explain the evolution of the concentration parameter with redshift. We also 
review some arguments used to derive the logarithmic slope of the inner and outer density 
profile. 
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1 INTRODUCTION 

Most of the mass in our universe is thought to be in the form of an 
unknown dark and collisionless matter that interacts with itself and 
with normal matter (baryons) only gravitationally. The formation of 
galaxies and the large scale structures in our universe is driven by 
the gravitational forces exerted by this dark matter (DM). On large 
scales the agreement of observations with models where the dark 
matter can cluster on all mass scales (cold dark matter: CDM) is 
very good. In this cosmological model, small-mass structures col- 
lapse first while larger galaxies and clusters form later by accretion 
and merging of smaller subclumps. In recent years, it has become 
possible to investigate the nature of DM in more detail both because 
of improvements in numerical simulations and because it has been 
realised that the Local Group of galaxies is a perfect observational 
laboratory in which to test the predictions of DM simulations. 

Observations of rotation curves in nearby dwarf galaxies and 
low surface brightness (LSB) galaxies generally show that their 
DM density profile has a flat inner core jvan den Bosch et al. 2000; 
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I de Blok et alj|200ll: IdTBlok & Bosmd 12002.: .Borriello & Salucci. 
l2001t ISalucci & Burkerd 120001) . This appears to be in disagree- 
ment with the results of numerical simulations of CDM haloes. 
Using A'^-body simulations, several groups find that halo density 
profiles generally have an inner cusp with slope a = 1 — 1.5 and, 
most remarkably, that the mean shape of the profile is universal: 
it is the same from dwarf galaxies to c lusters at any redshift (e.g. , 
iNavarro et alJl997tlMoore et alJl999l but see lSver& White! 19981: 



iKravtsov et ai]ll9 98: Subr amanian et al. 2000'; lina & Suto 2000). 
A large number of solutions to resolve the disagreement between 
observations and theory have been proposed. They can be divided 
into two main categories: (i) feedback effects of the baryons and 
stars on the DM density profile and (ii) modifications of the phys- 
ical properties {e.g. , scattering cross section, temperature) of th e 
DM particles (e.e.. lBode et'al]|200ltlSpergel & SteinhardJl2000h . 
The first class of solutions does not seem to be very effective in 
removing cusps and while the second has achieved moderate suc- 
cess in explaining the absence of a large population of dwarf haloes 
around the Milky Way, this comes at the expense of introducing ad- 
ditional inconsistencies with the standard cosmological model. 

The disagreement between observations and CDM simula- 
tions mainly concerns low-mass haloes and small length scales. In 
hierarchical scenarios, small mass haloes are mostly assembled at 
high-redshift. On the other hand A^'-body simulations have focused 
on studying the formation of Milky-Way type haloes or clusters 
(1O"M0 <M< 10" Mo) at low-redshift (0 < Zvir S. 3). Re- 
cently, I Ricotti 2003a, hereinafter R03) have used TV-body simula- 
tions to study the DM profiles of haloes with masses between 10* 
M© and 10^^ Mq, comparing them at their characteristic redshift 
of formation (0 ;S ^vir ^ 10). This is done using simulations with 
different box sizes and comparing the DM profiles of the haloes at 
different redshifts, when the most massive haloes in the simulations 
are composed of about the same number of DM particles. The re- 
sults show that the density profile of DM haloes is not universal, 
presenting cores in dwarf galaxies and steeper cusps in clusters. 
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This result could be a resolution of the conflict between the results 
of CDM simulations and the observed dwarf galaxy rotation curves, 
if we assume that the profile is stable and remains flat from 2: = 10 
(when the simulation is halted because of the small volume of the 
box) to low redshift. 

It should be mentioned that in a recent work, IColin et aU 
y003|) have analysed simulations with similar resolution and box 
sizes as those in R03 and find no evidence for a trend of the in- 
ner slope with the mass of the dark halo. This disagreement may 
be due to differences in their analysis of the simulations, as they 
have analysed the slopes of the haloes only at redshift z = 3 for 
all their simulations, irrespective of box size. It is well known that 
the smaller the box size the earlier a simulation ceases to be a rep- 
resentative volume of the universe. In a simulation with lh~^ Mpc 
(32 h^^ Mpc) box size at z ~ 10 (z = 3), modes with wave- 
length comparable to the size of the box become nonlinear and, for 
instan ce, the halo mass function becomes unreliable. Colin e t al. 
yOOJ) conclude that the results of their simulations are in disagree- 
ment with ROB because the simulations are identical; however, we 
note that they have analysed their haloes at redshift 2 = 3 while 
the haloes in R03 were analysed at redshift z — 10. As a conse- 
quence of the additional evolution from 2 = 10 to 2 = 3, their 
haloes are about ten times more massive. It is therefore likely that 
the disagreement is related to their attempt to analyse haloes at low 
redshift in a small cosmological volume. Indeed, R03 also finds that 
the cores observed in the ten most massive haloes at 2 < 10 even- 
tually evolve to cusps. Any possible trend of the dark halo inner 
slopes with halo mass can easily be confused with a mass depen- 
dence of the concentration parameter. In order to demonstrate that a 
trend exists, particular care must be taken to minimise any possible 
systematic errors, for example, by comparing haloes in simulations 
with different box sizes only when the most massive h aloes have 
compa rable number of particles (R03). More recently, ICen et all 
J2004|) have confirmed the results found by R03; in addition, they 
identify a redshift dependence of the typical halo profile. 

In addition to t he disagreement on the value of the inner slope 
jPowe^^lJ J2003h find q = 1 while a = 1.5 is advocated by 
iMoore et alj |l999)) other groups have also found deviations of 
the DM profile from a "universal shape". These deviations mani- 
fest themselves as dependences of the profil e shape on the mass , 
envirormient, cosmology or merging history j Sver&_WhitdU998 ; 
Kravtsov et ai]ll998l : ISubramanian et alJIZOOOl: Ijing & SutolboOC ; 
Jins 20005. In recent years, it has become possible to observe the 
DM profile of haloes at the two extremes of the halo mass function: 
in dwarf spheroidal galaxies with M ~ 10* M© and in clusters 
with M ~ 10^^ Mq. In §|2|we summarise observational evidence 
that appears to confirm the results found by R03, albeit with signif- 
icant uncertainty. 

The simulation results in R03 strongly suggest that the mass 
function of the accreting satellites determines the inner slope of 
the DM profile. However, the results might also be interpreted as 
a trend with dynamical state in the assembly process of haloes of 
different mass. In this paper we study the dynamical stability of 
the flat profiles of dwarf-size haloes at 2 = 10 by evolving them 
from their redshift of formation to lower redshift without further 
merging. We extract a few haloes at 2 = 10 and we evolve them in 
isolation using a direct TV-body integration code (see §|3|for details 
on the method). In §|4]we show that the profiles are dynamically 
stable and retain an almost flat core. As shown in §|5|the DM pro- 
file is consistent with the observed velocity dispersion of stars in 
the inner regions of the Draco dSph. Simple scaling arguments are 
given in §|6|to understand how the concentration parameter evolves 



with redshift and what determines the slope of the outer profile. We 
summarise the results in §0 



2 NON-UNIVERSAL DM PROFILE: OBSERVATIONS 

There is growing observational evidence that the DM density pro- 
file in the central regions of low-mass galaxy haloes is not univer- 
sal and that cored profiles are favoured over cuspy halo models. 
Despite the on-going debate about the correct interpretation of the 
gas rotation curves of low-mass disk galaxies, it is clear that these 
galaxies display a wide range of inner halo profiles. What is also ap- 
parent is that the data generally favour logarithmic density slopes 
close to 0.2 (e.g.. lde Blok et al] l2003V The gas-free Local Group 
dwarf spheroidal galaxies provide a cleaner test-bed for DM mod- 
els as they are often completely DM dominated at all radii and it 
is therefore possible to measure their DM content by treating their 
baryonic content as a massless tracer population. The recent iden- 
tification of a kinematically cold stellar substructure in the Ursa 
Minor dSph (K levna et alj 2003) strongly suggests that the dark 
halo of that galaxy has a central core. Using TV-body simulations 
iKlevna et alj y003) showed that this cold structure would survive 
for less than 1 Gyr if the DM core were cusped. Only if the DM core 
has a uniform density cor e can the cold sub structure survive for a 
Hubble time. Additionallv. lMagorriaiJ fcoOSh found a = O.SStoJa 
for the Draco dSph. 

Gravitational lensing observations sh ow that the inner sl ope of 
the DM profile in clusters is 1 < a < l.S. lDahle etai]^2003^ stud- 
ied six massive clusters with masses M ~ 10^* M© at 2 = 0.3. 
They find a = 1.4+°? and c = 1.7^5:3 ^^ the 68 % confidence 
level. ISmith et alj y00l|) finds a = 1.3 for the lensing cluster 
A383 that has a mass M ~ 10^* M© (2 = 0.188). The lens- 
ing cluster MS 2137-2353 has been studied independently by two 
groups. Their conclusions are quite different. Gavazzi et al. ( 2003|) 
find 0.7 ^ Of ^ 1.2, but they note that if the identification of a 
fifth lensed image is confirmed by HST observations, an isother- 
mal profile with a flat core fits the data better. ISandet alj 12003) 
find a best-fit a — 0.35 and a < 0.9 at 99 % confidence level 
for the same cluster. These results are exte nded and strength ened 
by the analysis of six clusters presented in lSandet alJ yOOJ). X- 
ray studies based on Chandra observations find that the inner slope 
of the DM profile is consistent with a — 1 i Arabadiis et al. 2OO4). 
lUewis et alJ '2003) use a generalised profile to fit their observations 
and find a f=s 1.2 ± 0.04 for the central slope of A2029. 



Thus there is some observational evidence for a dependence 
on halo mass of the inner slope of dark matter haloes, and hence 
for the non-universality of the dark matter profile, in agreement 
with the findings of R03. 



3 NUMERICAL METHOD 

We extract three haloes (HALO 1, 3 and 7) from the simulation 
with I/box = 1 h^^ Mpc in R03 at 2 = 10 and we evolve each of 
them as an isolated TV-body system using a direct TV-body code as 
explained in the next subsection. We can neglect the Hubble expan- 
sion since for r <C Vc{r) / H (z) ~ 10 kpc, where K is the halo 
circular velocity, the effect on the inner profile is negligible. The 
location of the dark matter core radius, rs ~ 200 pc, is not affected 
by the Hubble expansion. 
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3.1 Af-body simulations 

We follow the evolution of the low mass haloes using the direct A'^- 
body integration code NB0DY4 ( Aarseth 1999). In order to avoid 
spurious two-body interactions (that are non-physical for dark mat- 
ter particles) we have adopted a standard force smoothing at small 
particle separation with the same Plummer softening parameter, 
e « 11 pc (physical) as that used in R03. We have also run the same 
simulations without force smoothing to investigate the possible im- 
portance of spurious two-body effects, finding that the profiles at 
early times are very similar. Even at later times there are only neg- 
ligible differences between the evolved profiles in the two sets of 
runs (softened and un-softened). Initially, the three simulations per- 
formed contain, respectively, 54341, 46935 and 26565 particles of 
mass 7050 M0 and their evolution is followed for about 5 Gyr (12 
Gyr for the un-softened runs). The simulation s were performed on 
GRAPE-6 special-purpose computer boards JMakino et alJll997l) 
at the American Museum of Natural History, New York and at the 
Institute of Astronomy, Cambridge. At regular intervals during the 
integration the positions and velocities of all particles in the simula- 
tion were recorded and the density profiles calculated. The density 
profiles are calculated as in R03. The centre of the halo is calculated 
iteratively for particles within a given radius. The radius is reduced 
until the centres calculated using the centre of mass, density max- 
imum and the minimum of the gravitational potential, converge. 
The results converge within 2e, where e ~ 11 pc is the Plummer 
softening parameter. 



4 RESULTS 

The main aim of this paper is to investigate whether the profiles 
found by R03 are in equilibrium. In order to check the stability of 
the haloes, however, we do not need to evolve the halo for 12 Gyr 
(the time from z = 10 to z = 0) but for a much shorter time 
scale: namely the crossing time icr ~ 0.5 — 1 Gyr. On this time 
scale two-body relaxation can be safely ignored since icr < irciax, 
where irciax is the two-body relaxation time scale. 

Fig.Qshows the the circular velocity Vc = (GA/(r)/r)^ " of 
the thi-ee haloes with A/ ~ 2 - 4 x 10* M© virialised at 2: ^ 10 
and evolved in isolation for about 5 Gyr. The profiles are shown 
at the times t after z = 10 indicated by the labels. The vertical 
dashed lines show the location of the virial radius when the haloes 
formed at z — 10. Note that at z = the halo virial radius (defined 
as the radius where the halo mean density is about 200 times the 
IGM mean density) will be about ten times larger (see S I6.1> but a 
fraction ^ 20% of the halo mass will reside in the inner 1 kpc. If 
instead the halo is accreted by a larger galaxy the outer envelope of 
the halo will be lost and the mass and radius might be reduced to 
the values at the epoch of formation (z — 10). The profiles shown 
here are reliable in the range of radii r > 2e = 22 pc, where e is the 
force smoothing length and r ^ 1 kpc. At larger radii, the Hubble 
expansion and the mean density of the intergalactic medium (IGM) 
should be included to determine the evolution of the outer density 
profile. We will return to this issue in ii l6.II 

As the figures show, after about 0.1 Gyr, the profiles settle into 
a stable configuration and the central regions of the haloes (within 
about 1 kpc) do not evolve significantly during the next 5 Gyrs. 
The final density profiles of the three haloes are very similar to 
each other and have almost flat cores with inner slopes in the range 
Q ~ ~ 0.4. After f ~ 5 Gyrs two-body relaxation may start to 
be important in driving the evolution of a mildly cusped core into 
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Figure 1. The ciirular velocity Vc of 3 haloes with M ~ 4 X 10* M0 
virialised at 2 ~ 10 and evolved in isolation for a time t ^ tcr- The 
circular velocity profiles are shown at the times after z = 10 shown by the 
labels. Three haloes are shown: HALO 1 (M ft; 3.8 X 10* Mq), HALO 3 
(M Si 3.3 X 10* M0) and HALO 7 (M « 1.9 X 10* M©) from the 
^box = 1 h^^ Mpc simulation in R03. The dashed Hne shows the location 
of the virial radius when the haloes formed at z ^ 10. The lower limits on 
the X axis corresDond to twice the force softenine lensth. 
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iroiax = (2.6 Gyr) 



a hal o with a flatter core (see e. g.. iHavashi et al]|2003t Ijin et alJ 
l2003h . This can be observed in HALO 3 that aUt ^ I Gyrs has 
settled to an equilibrium slope with a ^ 0.4 and at t > 5 Gyrs has 
evolved to an almost flat core (a ~ 0). 

The two-body relaxation time is tioiax oc a^ /{m-pp), where a 
is the velocity dispersion of the particles, m^ is their mass and p is 
their density. Taking a ~ K and the core density for p we have, 

M y/' / 7050Mq \ / rco.e \ ''^ 
IQ»Mq) \ TUp ) ^200pcy 

Since both the density profile at i = and the velocity disper- 
sion o are approximately constant with radius, troiax does not de- 
pend strongly on the radius. For the haloes in our simulations the 
estimated two-body relaxation time is iroiax ~ 5 Gyr, consistent 
with the observed evolution of HALO 3. The equilibrium profiles 
of HALO 1 and HALO 7 (after t ~ 1 Gyrs) are almost flat and 
therefore the cores in these haloes do not arise from any two-body 
effects in the present work. Moreover, since the Hubble time at 
z — 10 is ten times shorter than trciax '^ 5 Gyr, two-body relax- 
ation effects are not important even in the R03 simulations from 
which the initial conditions where taken (see R03 for extensive dis- 
cussion of this issue). 

Even the later evolution (at t ^ 12 Gyrs) is not strongly af- 
fected by two-body relaxation, in the sense that the profile does 
not evolve significantly towards a post core-collapse cusp. Indeed 
one would expect that relaxation would turn our initially cored 
profiles into cusped profiles (via core collap se) on tim e scales 
t ~ 17trciax ^ 85 Gyr for our haloes (Taka hashilll995h . We do 
not see the effects of core collapse in our systems during the first 
12 Gyr even in the simulations that do not include force softening. 



5 COMPARISON WITH OBSERVATIONS 

In recent years, the volume of kinematic data available for Lo- 
cal G roup dSph galaxies has increased significantly (e.g.. lMated 
Il997r. Klevna et al. 2002. 2003) . These data make it possible for 
the first time to map the velocity dispersion profiles of these galax- 
ies and have confirmed that in at least one case, the Draco dSph, the 
high c entral velocity dispersion persists to large radii ( Klevna et al. 
I2OOII) . It is instructive to compare the observable dispersion pro- 
file expected for a dwarf galaxy whose halo mass distribution is 
matched by one of those simulated in the previous section. We as- 
sume that the projected light distribution in the dwarf follows a 
Sersic profile 



E,(i?) = Eo exp [- (i?/i?s 



(1) 



where R^ = 177 pc and ris — 1.2 which is a good fit to the in- 
ner surface brightness profile of Draco ( Klessen et al. 2003). If we 
further assume that the stellar velocity distribution is everywhere 
isotropic then integration of the Jeans equations yields the expected 
(projected) velocity dispersion as a function of radius 



a\R) 



/-fp.(^)K^x)(:r^-it:^)V^ 
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where ps is the density of the light. Fig.|2|compares the dispersion 
profile for HALO 1 at z = (s olid line) with the observed disper- 
sion profile for Draco JKlevnaetal. 2002) . The other lines show 
the dispersion profiles using synthetic scaled versions of HALO 1 
(that has a final mass of Afdm — 10* M0). In the top panel we 
show two haloes with the same core density as HALO 1 but with 
twice its mass (dashed line) and four times its mass (dotted line). 
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Figure 2. Comparison of the projected velocity dispersion profile for 
a Sersic light distribution embedded in the dark matter distribution of 
HALO 1 (solid line) with the observed profile for the Draco dSph (points 
with errorbars). In the top panel we also show two synthetic scaled versions 
of HALO 1 {Aldm — 10* M0) with the same core density but with twice 
its mass (dashed line) and four times its mass (dotted line). The dotted line 
in the bottom panel shows a halo eight times more massive than HALO 1 
but with h alf its core density (i.e., virialised at 2; < 10). Observed data 
taken from lKlevna et alJ J2002i) . 



The dotted line in the bottom panel shows a halo eight times more 
massive than HALO 1 but with half its core density. This profile 
provides a better fit to the overall shape of the dispersion profile. 
Haloes with total mass Mdm ~ 2 — 8 x 10* M^ provide reason- 
able fits to the observed profile. We note that lStoehr et alJ i2002h 
also found haloes in their simulations which were able to match 
Draco's dispersion profile. However, given that the resolution limit 
(in physical units) of those simulations was comparable to the edge 
of the light distribution in Draco it is not clear that their agreement 
with the observed dispersion profile is physically significant. 



6 SCALING ARGUMENTS 

In this section we present a short review of previous semi-analytic 
arguments in order to illustrate the physical mechanisms that de- 
termine the shapes of dark matter haloes. Using very simple argu- 
ments, we attempt to clarify the links between our numerical re- 
sults and theory. Using A^-body simulations R03 have shown that 
the density profile 



Pdm « 



X"{l + X)-'-'^ 



(3) 



1.3 + (M]/« 



with X = caR/Ra, a = (9 + 3n)/(5 + n) ~ ^.^ -,- v-Mm,i4 

l)/(iV/dm 14 + 1) 3nd constant concentration parameter ca ~ 7 
provides a good fit to all recently-virialised DM haloes from dwarf 
galaxies to clusters. Here, n is the effective spectral index of the 
initial power spectrum of density perturbations, -RA(A'fdm) is the 
virial radius and Afdm,i4 = Afdm/(3 x lO^'* M0), where A/dm is 
the halo mass. 

The relationship between the inner slope a and n can be un- 
derstood using simple scaling arguments (see S I6.I> . In the present 
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paper, the haloes are evolved from z = 10 to low redshift in a 
zero-density universe, neglecting the Hubble expansion. This sim- 
plification prevents us from knowing the evolution with redshift of 
the concentration parameter and the outer profile of the dark halo. 
In S I6.2l we review some simple arguments whose aim is to clarify 
the physical processes that determine the evolution of the concen- 
tration parameter and the outer slope of the density profile from the 
time of halo virialisation to 2; = 0. 



6.1 Inner slope 

In this section, following IPeebleJ il974) and lHoffman & ShahamI 
il985l) . we review the arguments used to derive the inner profile 
of DM haloes under the simplistic assumption that each sub-halo 
(here represented as a step function in density) survives undisrupted 
when incorporated in a larger halo. The result that we obtain agrees 
remarkably well with the simulations in R03. 

Assuming an Einstein de-Sitter universe, linear perturbations 
5k grow as 5k{t) oc a{t) (x f^'"^, where a is the scale factor of 
the universe and t is the age of the Universe. The variance a^ in 
the amplitude of the density perturbations varies with perturbation 
mass M and time t as (j'\M,t) ex f*/3Af-(=^+")/3. Here n is 
the effective spectral index of the initial power spectrum of density 
perturbations P{k) oc fc". When cr^(M, i) reaches approximately 
unity, the associated perturbations become non-linear: this occurs 
at time fvir(M). Hence, tvir oc M^"^^'"^: small masses become 
non-linear at earlier times if n > ~3. The halo sizes as a function 
of mass are given by R{M) = a(tvir)rco oc M^"^'"''^^, where 
the comoving radius is r^o oc M^'^ , and the mean densities of 
the haloes are p{M) oc a(tvir)"^ oc M"(^+"^/^. The profile of 
a halo formed by the merger of many such haloes is then given by 
the envelope that contains all the merged haloes piled on each other 
like "Russian dolls". Smaller, higher-density haloes will reside in 
the centre while the more extended, lower-density haloes dominate 
the profile at large radii. Combining the expressions for R{M) and 
p{M) we find p{R) oc R°', where a = {9 + 3n)/(5 + n). The 
effective n in the 1 h^^ Mpc simulations of R03 is —2.6 ± 0.1 
corresponding to a = 0.4 — 0.6. 

6.2 Halo concentration and outer slope 

The post- virialisation evolution of isolated dark haloes can be un- 
derstood (to a first approximation) using simple argu ments partially 
based on t he theory of secondary infall developed bv lOunn & Gotll 
mzlandlGotJIliTl. 

Particles that decouple from the Hubble flow and become 
gravitationally bound are found, by definition, at radii where the 
escape velocity is smaller than the Hubble expansion (i.e., Vcsc ^ 
'^Hubble = H{z)Riis_io, where H{z) is the Hubble parameter). 
Therefore, since Uesc ~ K, where K is the circular velocity, we 
have 



K = 



GM 



HRhalo, 



(4) 



-Rhalo 

where i?haio is the outer radius of the dark halo. Since the Hub- 
ble expansion decreases monotonically with time and the mass of 
the halo after formation is either constant or increases with time, 
it follows that a halo that evolves in isolation (i.e., without major 
mergers) increases its size and mass. Assuming, for simplicity, that 
the halo mass increases with time as M{t) oc t^, from equation 
we have _Rhaio oc ii'~^/^M(f)^/^ oc //~(^+^'''^ or, neglecting the 
cosmological constant, _Rhaio (2) /-Ra = [(l + 2vir)/(l + 2)]''"^^ , 




Figure 3. (top) Slope a of the inner DM profile as a function of halo mass, 
A^dm (solid curve). The points with error bars are based on observations 
(see §|2|. (bottom) Concentration parameter, ca, of the generalised profile 
at z = (solid line), z = 3 (dashed Une) and 2 = 10 (dotted line) as a 
function of halo mass. The points with errorbars are the concentrations at 
z = for the same haloes as in the top panel. We estimate the virial radius 
of dSphs assuming that the DM density decHnes as p oc R^'^ outside the 
stellar tidal radius. Note that Sand et al. (2004) find an almost flat slope of 
the inner density profile in clusters, in disagreement with Dahle et al. (2003) 
and others. 



where Ra is the virial radius at the time of formation. According to 
the numerical results presented in §|4|the dark halo core radius, rs, 
remains constant with time. Therefore, the concentration parame- 
ter, CA = i?haio/rs, evolves as 

CA{z,Z^ir) = CA(2vir) ( , , ^"^ ] ■ (5) 



1+2 J 

Here, we note the result found by R03 that the concentration pa- 
rameter CA(2vir) of recently virialised haloes is roughly constant 
and independent of the halo mass. Finally, based on simple scal- 
ing arguments (see GotI 1975), the outer envelope should have a 
dens ity profile p oc M(2)iihaio(z)"^ <x R-^^^^+^\ 

iGunn & Got3il9721 find that during secondary infall ^ = 2/3 
(infalling shells do not cross during the infall). But if, after the col- 
lapse, crossing shells start leading to violent relaxation then one 
expects the equilibrium halo to be more centrally condensed and 
the outer envelope should have a logarithmic slope between —4 
(violent relaxation) and —2.25 (infall upper limit). We also note 
that the assumption that the infall rate is a power law in time is an 
oversimplification since isolated haloes will evolve differently de- 
pending on whether they are located in an overdense or underdense 
region of the IGM or as a consequence of extended density cor- 
relations. In other words, the density power spectrum introduces 
preferred length scales, and local departures from fim = 1 in- 
troduce preferred time scales. This more complicated and realis- 
tic accretion history is reflected in a large scatter of the concentra- 
tion parameters o f haloes of a given mass, as shown by N -body 
simulations (e.g.. iBullock et alj|200il: IWechsler et alj|2002h . For 
the purpose of this section we simply note that assuming ^ ^ 
(i.e., neglecting the mass accreted after virialisation) we obtain 
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good agreement with TV-body simulations for the dependence of 
the concentration parameter on the redshift and formation redshift, 
ca oc (1 + Zvir)/( 1 + z) (e.;^.,[W echsler et al. 2002) and outer pro- 
file p oa r~^ (e.g., [Navarro et al. 1997). Our simpl e arguments are 
also in good agreem ent with the numerical results of lWechsler et alj 
J2002hnzhaoetalj t2003bia) that find that the concentration of 
haloes is about constant (cs ~ 4) at the end of an initial phase 
of rapid accretion (i.e., virialisation) and increases later to larger 
values as a consequence of a slower secondary phase of accretion 
(i.e., secondary infall). In agreement with those numerical studies, 
equation J5j implies that at high redshift most haloes have small 
concentrations ca ~ 3 — 4, since they did only recently virialise. 

In CDM cosmologies small mass haloes form earlier. We can 
use the results of S I6.1l to relate the mass of the halo to its typ- 



ical redshift of virialisation: (1 + Zvir) oc M^ 



■(n + 3)/6 



Using 



the relation a — {9 + 3n)/(5 + n) and ca oc (1 + Zvir) we 
have ca{z = 0, M) oc M^^ ~ "' . In Fig . Owe compare ca 
and Q as a function of mass with observations JKlevna et all2003l : 
iMagorriaijEoOSI : Ide Blok et alJl2003l : iDahle et alJl2003h . The DM 
profile slope a and concentration ca based on the generalised DM 
profile are consistent with the observational data if the concentra- 
tion parameter at virialisation is CA(zvir) ~ 3 — 4 ^. Note that the 
concentrations shown in the bottom panel are a rough approxima- 
tion because they have been derived using scaling arguments and 
are valid only qualitatively. At high redshifts the halo concentra- 
tions are about constant and almost independent on the halo masses, 
as illustrated by the dashed (at z — 3) and dotted (at z — 10) lines 
in the bottom panel of Fig.|3| 

Finally we note that the mass-dependent profile shape coupled 
with a mass-independent concentration parameter at virialisation 
gives first-order results in good agreement with a universal profile 
shape coupled w ith the m ass-dependent concentration parameter, 
cnfw, found by JBullock et al. 2001) and others. This renders the 
results presented here more consistent with other results than one 
might first have guessed (for a discussion see R03). 



of the dark haloes in our simulations resembles the observed profile 
in the inner regions of Draco dSph, making these haloes plausible 
hosts for the observed dSph galaxies. 

In this paper we have evolved the simulated haloes as isolated 
systems. While some of the satellites of the Milky Way (e.g., the 
Sagittarius dSph) are clearly in the process of being destroyed by 
the tidal field of the Milky Way(e.g., Mateo 1998), others (e.g., the 
Draco dSph) show l ittle evidence for signifi cant tidal ef fects in their 
inner regions (e. e. . iKlevna et alj.2002: .Klessen et aljt2003.) . It is 
therefore reasonable to assume that their inner haloes have evolved 
in near-isolation since the formation of their stellar populations. 
The persistence of the central DM core in our simulations confirms 
the conclusion of R03 that the mass dependence of the inner slope 
of dark matter haloes could arise solely from the power spectrum 
of primordial density fluctuations. 

The concentrations and slopes of the inner profiles of sim- 
ulated haloes are found to be in good agreement with avail- 
able observations of dwarf galaxies and clusters at 2 = 0. The 
concentrations as a function of mass and redshift agree with 
those found in published A'^-body simulations in the mass range 
10^^ M0 ^ M ^ 10^" Mq covered by those simulations. We have 
also summarised scaling arguments which demonstrate why the in- 
ner slope depends on the halo mass and why the concentration pa- 
rameter depends on the time elapsed since virialisation. 
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7 SUMMARY AND CONCLUSIONS 

Cosmological simulations of the formation of the first galaxies at 
z > 10 show a good agre ement with the observ e d properties of 
dSphs in the Local Group 1 Ricotti et al. 2002a b; Ricotti "2003b|; 
iTassis et al. 2003 : Susa & Umemura 2004) . Those simulations, be- 
cause of the inclusion of radiative feedback effects, can repro- 
duce the observed number of luminous galactic satellites and the 
observed internal properties of dSphs (e.g., mass-luminosity rela- 
tion, metallicity-luminosity relation, mass-to-light ratio, half-light 
radius). Acc ording to those simulations and similar pure A'^-body 
simulations iRicottill2003ah . the DM profiles of dwarf galaxies at 
z = 10 have a central core instead of a steep cusp. 

In this paper we have used A'^-body simulations to show that 
the central DM core found in dwarf galaxy haloes at z — 10 is 
dynamically stable. At low redshift the halo still has a flat core 
in agreement with the observed rotation curves of low-mass disk 
galaxies and evidence from the Ursa Minor dSph galaxy. We find 
that the projected velocity dispersion profile for a dwarf galaxy with 
a light distribution similar to that of the Draco dSph residing in any 



^ Note that the haloes in R03 have ca ~ 6— 7, but they have been analysed 
some time after virialisation. A 40 % increase of their concentration since 
virialisation is plausible. 
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